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Use of a cDNA microarray to 
analyse gene expression 
patterns in human cancer 

Joseph DcRisi 1 *, Lolita Penland 2 & 
Patrick O. Brown 2 (Group 1); 
Michael L. Bittner 3 *, Paul S. Mehzer\ 
Michael Ray 3 , Yidong Chen\ Yan A. Su 3 & 
Jeffrey M. Trent 3 (Group 2) 

The development and progression of cancer 1 " 3 and 
the experimental reversal of tumorigenicity 4 - 5 are 
accompanied by complex changes in patterns of 
gene expression. Microarrays of cDNA provide a 
powerful tool for studying these complex phenom- 
ena 6 " 8 . The tumorigenic properties of a human 
melanoma cell line, UACC-903, can be suppressed 
by introduction of a normal human chromosome 6, 
resulting in a reduction of growth rate, restoration 
of contact inhibition, and suppression of both soft 
agar clonogenicity and tumorigenicity in nude 
mice 4,5 - 9 . We used a high density microarray of 
1,161 DNA elements to search for differences in 
gene expression associated with tumour suppres- 
sion in this system. Fluorescent probes for 
hybridization were derived from two sources of cel- 
lular mRNA {UACC-903 and UACC-903(+6)] which 
were labelled with different fiuors to provide a direct 
and internally controlled comparison of the mRNA 
levels corresponding to each arrayed gene. The flu- 
orescence signals representing hybridization to 
each arrayed gene were analysed to determine the 
relative abundance in the two samples of mRNAs 
corresponding to each gene. Previously unrecog- 
nized alterations in the expression of specific genes 
provide leads for further investigation of the genet- 
ic basis of the tumorigenic phenotype of these cells. 

DNA microarrays. containing Lift I total elements, 
including S70 different cONAsand controls*" 11 (see 
Methods), were printed robot ically onto a glass micro- 
scope slide in four quadrants covering an area of about 
1 cm 2 (Fig. I ). We prepared fluorescent cDNA probes 
using total poly ( A) + mRNA from UACC-903 cells and 
UACC-903(+6) cells by labelling with a green and red 
fluor, respectively. A mixture of the two flourescently 
labelled probes was hybridized to the DNA microarray. 
This comparative hybridization method, coupled with 
the doping of synthetic standards and an estimation of 
statistically significant deviation for local background 
variance allowed a direct and quantitative comparison 
of the relative abundance of individual DNA sequences 
in this complex sample'*"*. We added a set of synthetic 
poly (A) + -tailed 'mRNAs' to the purified mRNA from 
each cell line as internal standards to assist in quantita- 
tion and estimation of experimental variation intro- 
duced during labelling and reading. Targets 
complementary to these standards were included, in 
duplicate, on the microarray. Based on these standards. 
mRNA species comprising 1:10.000 of the mass of the 
polv (A)* RNA could readilv be delected. 

in a representative two-colour fluorescent scan of all 
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spond to genes preferentially expressed in the tumori- 
genic UACC-903 cell line, and the reddish spots corre- 
spond to genes preferentially expressed in the 
non-tumorigenic UACC-903(+6) cell line. Genes 
expressed at approximately equal levels in the two cell 
lines appear yellow or brown. A portion of the array at 
higher magnification highlights the diverse pattern of 
differential expression observed (Fig. 2b). In Fig. 2c, rec- 
tangles corresponding to specific array elements are 
coloured to reproduce the hue and intensity of the fluo- 
rescent signal at each element. The hybridization signals 
from a duplicated set of genes arc shown juxtaposed, to 
illustrate the reproducibility of the hybridization signals 
for each gene. 

To address the possibility that an apparent difference 
in expression might result from experimental variables 
unrelated to the difference in chromosomal composi- 
tion between the two cell lines, we examined the vari- 
ance in expression for 90 'housekeeping* genes. We 
selected these genes based on the assumption that they 
would not be differentially expressed between the two 
cell lines. The averaged red/green ratio for this subset of 
genes was 1.13. The averaged red/green ratio for the set 
of five internal standards was 0.97 (// = 10). The vari- 
ability in the expression level of the housekeeping genes 
probably overestimates the experimental variability in 
measuring differential expression. As a conservative stan- 
dard, an absolute fluorescent signal (red or green) with 
an intensity greater than that observed at the control 
array elements containing total human genomic DNA 
was considered to represent specific hybridization. Gene- 
specific hybridization was therefore only considered sig- 
nificantly different between samples if the following two 
criteria were met: i) the signal intensity (green or red) 
exceeded this threshold; and ii) the louarithm of the 
red/ureen fluorescence siunal ratio differed bv >3 S.D. 
from the mean logarithm of this ratio for the 'house- 
keeping* gene pane! (that is, ratios <0.52 or >2.4). 

By the.se criteria. mRNA levels for 1 5/870 ( 1 .7%) genes 
were significantly diminished, while the mRNA levels 
for 63/S70 (7.3%) genes were significantly increased in 
association with suppression of tumorigenicity by intro- 
duction of chromosome 6. To test the reliability of 
microarray hybridization results in identifying differen- 
tially expressed genes, we analysed 16 genes by north- 
ern analysis. In each case, the results of northern analysis 
corroborated the differential gene expression identified 
by microarray hybridization (Fig. 3). 

Significant differences in expression between these 
two cell lines identified several genes as candidates for 
determining features of the tumorigenic phenotype of 
the melanoma cells. For example, among the genes 
delected with significantly higher expression (> 10-fold) 
in the tumorigenic cells was the human brown locus pro- 
tein (TRIM/melanoma antigen gp75). This is the most 
abundant glycoprotein in melanocyte cells and a critical 
melanosome membrane protein 12 * 1 -*. Additi nally, its 
expression is reduced when melanoma cell lines are 
induced to differentiate by treatment with HMDA 12,13 . 
Also expressed at a significantly higher level was a spliced 
variant of the mRNA encoding myelin IM .IVPM20. This 
is widely expressed in neural crest derived cells in early 
development and has been suggested to play a role in 
cell-cell signaling during development 1,1 . 
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• (54) Hybridization Specificity Controls 
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Fig. 1 Properties crt cONA microarrays. *. A fluorescent scan of DNA printed onto a poJy-rysme coated slide. The DNA is stained with a DNA-specHjc fluorescent 
dye. YOYO. The center- to-center spacing of adjacent spots is 450 u. allowing the potential tor up to 10.000 spots/2.54 X 7.62 cm microscope slide. 6. Effi- 
cient blocking of hybridization to DNA repeats. Hybridization of fluorescein-iabetled poly (dT)* to arrays in the absence of competitor produces strong 
hybridization to immobilized poly (dA)* as well as to some cDNAs. such as the EST T64827 shown. Rhodamine- labelled cDNA (red) from the UACC-903 cell 
line hybridized in the presence of poly (dA)' blocker shows little if any signal at either srte (Total H = total human). Similarly, hybridization with tiuorescetn-iabetted 
Con DNA in the absence of competitor produces bright signal on immobilized Cotl DNA. total human DNA and at some cDNA elements (presumed to con- 
tain highly repeated sequences, such as R23416): while Rhodamine-labetied cDNA (red) from the UACC-903 cell line produces little if any signal at these 
locations when hybridized m the presence of excess unlabetled poly (dA)'. and human Cot 1 DNA. The absence of signal at some cDNA locations following 
UACC-903 cDNA hybridizations also indicates that the PCR-amplified. plasmid vector sequences at ail cDNA targets oo not contribute significant hybridiza- 
tion signal, c. Schematic of the array organisation. Robotic printing from 96 well microliter trays was carried out with 4 print heads, spaced to fit into 4 adja- 
cent microliter wells. This maps the contents of each tray into four separate quadrants on the glass slide. A cotour-coded map of the general distribution of 
target types in each of the resulting quadrants is shown. 



els were devilled bv ihe addition of a normal chromo- 
some 6 ( 1 7 genes) are known to he activated by I FN -v. a 
cardinal proinflammatory cytokine that, among other 
activities, induces expression of the gene products of the 
MHC class II locus. For example, the mRNA encoding 
monocyte chemotaetic protein I (MCAF/MCPI ), a 
cytokine that induces monocyte chemotaxis and activa- 
tion ,, * ,,# \ was more than 1 0-fold less abundant in the 
tumorigenic cell line. In the skin, MPCl is critical in the 
regulation of cutaneous monocyte trafficking 1 *" 1 ', and 
elevated expression plays a role in suppression of tumour 
growth and metastasis 1 **"" 1 . The mechanism by which 
these interferon- Y regulated genes are induced in UACC- 
903 cells by transfer of a normal chromosome 6 remains 
to be determined. It is worth noting, however, that the 
interferon-? receptor gene is localized to the distal long 
arm of human chromosome 6. 

Finally, several genes that showed > 10-fold higher 
expression in the suppressed UACC-903(+6) cells have 
previously been recognized in other models of tumour 
suppression. Most notably, there was elevated expres- 
sion of the mRNA encoding WAF1 (p21 ), a key media- 
tor of tumour suppression by p53 (ref. 18). The p21 
protein had previously been identified as a melan ma 
difTerentiati n-ass ciated antigen (termed mda-6) ,9:o . 
In melan ma cell lines suppressed for metastasis by the 
introduction of chromosome 6, expression of WAF1 
< p2 1 ) mRNA and protein c rrelates inversely with 
metastatic potential* 0 . 



These results provide a wide view of the diverse sys- 
tems that are altered in this model system of tumori- 
gen icily, and focus attention on specific gene products 
and pathways that may be of particular importance in 
this tumour type. 

Our ability to classify human cancers in a way that 
reflects the underlying molecular pathology or that 
anticipates their potential for progression or response 
to treatment, remains primitive. Using cDNA microar- 
rays to define alterations in gene expression associated 
with a specific cancer may be an efficient way to uncov- 
er clues to the specific molecular derangements that con- 
tribute to its pathogenesis and thus identify potential 
targets for therapeutic intervention. Moreover, recogni- 
tion of pathognomonic alterations in gene expression 
might provide a basis for improved diagnosis and mol- 
ecular classification of cancers and thus allow selection of 
the most appropriate therapeutic strategies. 

Public databases of human expressed gene sequences 
contain partial sequences of at least 40,000 different 
human genes 11 , and efforts to develop a human tran- 
script map have developed rapidly- 1 . Based on the high 
yield of information obtained, using an array of < 1,000 
different genes, a more comprehensive survey of gene 
expression patterns, using a more complete array of 
human genes, will likely pr vide a rich source of new 
and useful insights into human biology and a deeper 
understanding of the gene pathways involved in the 
pathogenesis of cancer and other diseases. 



Fig. 2 DNA microarray analysis of changes in gene expression between the tumongenic cell line. UACC-903. and its non-tumorigenic derivative UACC903K6I 
derived by introduction of a normal chromosome 6. a. A ratio image of the results of simultaneous hybridization o! Rhoaamine110Mabelled cDNA loreeni 
from UACC-903 and Cy3-labe»ed cDNA (orange-red) from UACC-903C6) to a microarray. To produce this image, the scan ,mages corresponding to ^flu- 
orescent probe were combined as the appropnate colour channels in a single image. Arrows indicate the location w.tmn the array of the corresoondina 
genes analysed by northern blotting (Fig. 3). b. A magnified image of the area of the array boxed in white in (a), c. Maon.f.ed imaqe of three cDNAs indenti- 
hed by arrows ,nja). represent.ng the cDNAs for: left. MCARMCP- 1 (r/g rat.o > 1 0): centre, fractm <r/g rat.o 1 .04): and" right. »- 1 ^nuchymotrypssn (r/g ratio 
0.2) {see Fig. 3j. d. simphf^d representation of rat© hybridization results. Quantitative fluorescence intensity data is extracted from each array target "Die aver- 
age target colour ratio determines the hue of each box and the average intensity determ.nes the bnghtness of each box. in this image the ordeTof'the boxes 
corresponds to the.r original order in the m.crotiter plate from wh.ch they were printed. Duplicate pr.nt.ngs of the same plate can be examined side bv side 
as in the first two rows shown here, to assess reproducibility of the hybnd.zat.on results (see text) Numbered arrows indicate the location within the arrav cor- 
respond^ to genes analysed by northern blotting in Fig. 3. y 



Methods 

Generation of micruarrnvs. hybridization, scanning. The 
preparation ot mated micro**. ope .slides ami subsequent robot- 
it printing ol DNA was carried out in .i manna similar in thai 
described . Briefly, pre -v. leaned glas» slides were treated with 
poly. I. -lysine soliitmn t Sigma i to form ;in adhesive mi. lace f«>r 
printing. PC.K products, purified by ethanol purification, were 
resuspended in .V\ SSC A custom built arraying robot picked 
up and deposited small volumes {-5 nanohiers) of DNA onto 
the slides. After printing, the slides were washed in a SI >S 
sol til ion. The remaining bound DNA was denatured by sulv 
merging the slides in 95 "('.distilled water lor 2 min followed bv 
.i hriel wash with 95% ethanol. DNA was I'V c.osxlinked m the 
slides (Stratagcne Siratalinker. Wl ml). Tn prevent non-specific 
probe binding, the slides were blocked by rinsing in a solution 
of 70 mM succinic anhydride dissolved in 0.1 M boric acid pH 
K.n. containing 35% I -methyl -2- pyrrol idi none (Aldrich). 
Additional protocols and parts list pertaining to microarrav 
labrication can be obtained from http://cmgm.stanfonl.edu/ 
pbrown. 

Purified; labelled cDNA was resuspended in 1 1 pi of J.5x SSC 
containing 4 pg ol poly UlA)* DNA. 2.5 pg t.. ro/i tit N A. A pg of 
human Cot I I >NA ((iibco IIKI ). and II..* pi ol HIV.. SDS. Prior to 
hybridization, the solution was boiled for 2 min then allowed in 
cool to loom temperature. Hybridization wa> carried out ai 
rO Q C for - N h in a water bath. Prior to wanning, slides wete 
washed in > SM !. 0.:» M SI >S lor 5 min ami O.Jy SS< '. for I min. 

local laser microscope built bv S Smith with software written 



by N. V.w. A separate scan, using the appropriate excitation line, 
was done lor each ol the luo lluorophorcs used. Data was col- 
lecied .n a maximum resolution of 9 microns/pixel with 12 bits 
ol depth. 

Probe preparation and labelling. KNA was extracted from cells 
using thcTi ia/ol reagent I I 'l l Inc. I. following the manufactur- 
ers directions. cDNA probes were synthesized from singly oligo 
d I -selected ( Pharmacia i mKNA pools. Ruorcsccntly labelled 
el IN A was prepared from mKNA by oligo dT- primed polymer- 
ization using SuperScr.pt II reverse transcriptase UTI Inc.). 
The pool of nucleotides in the labelling reaction was 0.5 mM 
dcVIP. dATP and dCI'P and H2 mM dTTP. Fluorescent 
nucleotides. Khodamine III) dUTP (Pcrkin Timer Cetus) or 
CyJ dlTP lAmersham), were present at 0.1 mM. Probes were 
purified by gel chromatography (BinSpin h/Binllad) and 
ethanol precipitation. 

Selection of cDNA elements and generation of control tem- 
plates. Synthetic cDNAs were prepared by cloning random 
/inwHI and Hindi 1 1 ended fragments of E. toli DNA in the vec- 
tor pSPM poly (A)' (Promega). linearizing isolated plasmid 
DNA with fcmUi and synthesizing poly I A)' tailed KNA com- 
plementary to the insert l.om (he resident SIN* promoter 
(Piomega). Prior to use. the synthesized UNAs were selected on 
oligo d T cellulose. The largesi group ol cDNAs consisted of 674 
cDNA clones from the 1MB arraved normalized infant brain 
'■■■•a. i liiese . .v..„*icJ i«* ....»». •••cry NIB 

lihran member that corresponded to a named gene according 
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Fig. Z Nortnern nyor idnaiton suostamtai* 
mg the consistency of the cDNA mrcr oar- 
ray results. Corresponding locations wrfrten 
the cDNA microarray illustrated m Fig. 2a 
are provided tor 1) Waf-Vp2V. 2) MARCKS; 
3) collagenase: 4) MCAF/MCP-1; 5) rz-T- 
antichymotrypsin: and 6) 0 -acf/n. The sig- 
nal detected by a radio-labelled 0-acrm 
probe represents a control for loading vari- 
ance, with a red/green ratio observed on 
the cDNA microarray (Fig. 2a.c) for p-actin 
of 1.04. 

to the UniGene EST clustering sys- 
tem- 1 - : . The second kirpest itroup of 
clones con nisi eil of' 1 83 sequenced cDNA 
clones generated by subtraction of cDNA 
from I he chromoMimc-A suppressed 
non-iumorigenic UACC-903 (♦*) cell 
line with cDNA from its parental tumori- 
genie cell line UACC-903 ( ret. 9). 
Approximately MM) additional genes 
(total 870 genes arrayed) were obtained 
from KST libraries on the basis of their 
n- 1 -Antichymotrypsin expression pattern 1 1 issue specific, and so 

on I. hach array included the following 
hybridization controls: plasm id vector, 
lambda. OX 174 phage, total human PNA. 
human Coll UNA. and poly (A)\ The 
svnthetic standards used for normaliza- 
lion of signals in each wavelength were 
also arrived, ("out nils were included in 
each quadrant of the array to assess the reproducibility of the 
hvhridi/aiion signal. Two plates of cDNA clones (derived from 
the UAliOJOj subtracted library) were also arrayed in dupli- 
cate, fidelity of the Unigene array relative to dhKST was tested 
hy sequencing of a random sample of 1 1 clones used for 
microarray construction. All sequences were identical with the 



MCAF/MCP-1 



P-Actin 



corresponding dhLST entries. Additionally each mM.un.nvJ 
cDNA from the 1-ACC-W3 subtracted library wa< H-queiued. \ 
listing of cPNAs comprising this microarray w\u<\\ ueu 
derived from the Unigene and 'housekeeping panel »..m iv 
obtained from hltp://www.nih.gov/|)lU/LCr./AUKAVe\pn . 
html. 

Northern blot analysis. Total RNA, 10 pg per lane, was elec- 
trophorescd in 1.2% agarose* formaldehyde gels and transferred 
onto nylon membrane (Hyhond-N*. Amersham) by capillary 
blotting overnight. For DNA probes insert fragments hum the 
Scares 1Mb cDNA library 1 " were obtained by vector PCR for 
p2l, MAKCKS. (i* I -antichy nun rypsin anj act in. Htobcs lor 
fibroblast collagenase and NH*AF/\U!P-I were isolated from a 
UACC-Wtfl+fO enriched cDNA library" with all probe* 
labelled by random priming. Filters were washed to a stun- 
prncy of O.lx SSC at 42 *C for 20 min. 

Web sites. hitp:/ycmgm. stanford.edu/pbnnvn for protocols and 
parts list pertaining to microarray fabrication. 
htip://www.nchgr.nih.gov/l )IH/L(".Ci/AKKAY/expn.himl for a 
listing of cDNAs comprising this microarray which were 
deriml from the Unigene anil 'housekeeping" panel. 
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